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Abstract: Phenols as a class of molecules have been reported to exhibit higher log 
maximum fluxes through human stratum corneum, SC, from water, log /mhaq, than other 
classes of molecules. This suggests that their corresponding log maximum fluxes through 
silicone from water, log 7mpaq, may be useful to extend the existing n = 63 log 7mpaq 
database to include more log 7mpaq values greater than 0.0. The log 7mpaq values for n = 1 
phenols predicted to give log 7mpaq values greater than 0.0 based on their log /mhaq values 
have been experimentally determined. These n = 1 new log /mpaq values have been added 
to the existing n = 63 log 7mpaq database to give a new n = 70 database and the n = 1 
literature log ./mhaq values have been added to the existing n = 48 log /mhaq database 
(matched to the n = 63 log 7mpaq database) to give a new n = 55 database. The addition of 
the n = 1 phenols improved the correlations of these flux databases when fitted to the 
Roberts-Sloan equation, RS, as well as the correlation between the matched experimental 
(Exp.) log 7 M paq with the Exp. log 7 M haq- 

Keywords: maximum flux; silicone membrane; Roberts-Sloan equation; water solubility; 
lipid solubility; human skin in vitro 



1. Introduction 

The rate-limiting barrier to diffusion of molecules through human skin is the stratum corneum, SC. 
The SC is comprised of highly dense, polar, proteinous corneocytes embedded in a lipid matrix. The 
lipid matrix in turn is comprised of multiple lipid bilayers containing mostly ceramides, fatty acids and 
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cholesterol [1-3]. A tortuous pathway through this lipid matrix and around the corneocytes is generally 
considered to be the route followed by molecules diffusing through the SC. Since the route of 
molecules diffusing through the SC is lipid-like, the lipid solubility of molecules diffusing through the 
SC is an important physiochemical determinant of the efficiency of the diffusion process. Similarly, 
any surrogate for the SC must present with substantially lipid-like properties. Since silicone 
membranes present with highly lipid-like properties, it has been suggested that silicone membranes 
could be a surrogate for skin in diffusion cell studies and that those results could be used to predict 
diffusion through human skin [4,5]. Also, from a theoretical basis it can be assumed that, if the 
diffusion of unit mass per unit area per unit time, flux, through silicone membranes can be accurately 
modeled by an equation (Roberts-Sloan, RS, [6]) that accurately models maximum flux through 
human skin [7] exists, a linear relationship exists between the maximum flux of molecules through 
silicone membranes from water, log 7mpaq, and their maximum flux through human skin in vitro from 
water, log 7mhaq [5]. Thus experimental (Exp.) log 7mpaq could be used to predict Exp. log 7mhaq- 

The flux data for molecules from which the Roberts-Sloan, RS, equation was derived is based on 
the application of saturated solutions (suspensions of molecules in a solvent) to the membrane being 
used to give maximum flux, J M - Thus all the molecules are presented to the membrane at their 
maximum thermodynamic activity in that solvent [6-8]. Since the molecules are presented to the 
membrane at their maximum thermodynamic activity, at equilibrium the molecules are also at their 
maximum thermodynamic activity in the membrane [8,9], i.e., at their solubility limit or saturated 
solubility in the first few layers of the membrane, S M i- Thus, J M depends only on the solubility of the 
molecules in the first layers of the membrane, and unless the solvent changes the solubility limit of the 
membrane, the solvent/vehicle has no effect on J M . The form of the RS equation for predicting J M 
derives from an expansion of Fick's law, Equation 1, so that the dependent variables are molecular 
weight, MW, solubility in a lipid, Slipid, and solubility in water, S A q- Smi can then be estimated from 
the product of the partition coefficient between the vehicle (water in this case, AQ) and a surrogate 
lipid for the membrane (octanol in this case, Soct), (Koct.aq) 3 -constant, and the solubility in the 
vehicle, Saq: (£ocT:AQ) >- constant-S A Q. Expansion of that product into solubilities and taking the logs gives: 
log Smi = y log Soct - y log S AQ + log S AQ + log constant = y log Soct + (1 - v) log S AQ + log constant. 

J = D/L (C M1 — C Mn ):/ M = D/L (5 M1 — C Mn ) (1) 

where D is the diffusion coefficient of the molecule in the membrane, L is the thickness of the 
membrane, Cm\ is the concentration of the molecule in the first few layers of the membrane and Cmh is 
the concentration in the last few layers of the membrane which is assumed to approach zero. A linear 
relationship must exist between log D + log Sm\ of molecules in a silicone membrane and log D + log Smi 
of molecules in human skin in order for a linear relationship between log 7mpaq and log /mhaq to exist. 

One problem with determining if the linear relationship between log 7mpaq and log /mhaq exists is 
that there are only about n = 63 molecules for which log 7mpaq (output) and the necessary 
corresponding physicochemical properties (log Saq and log Soct, input) literature values exist which 
can be fitted to RS [5]. Of those n = 63, only 18 molecules exhibit high output values; log 7mpaq values 
greater than 0.0. Simple phenols present an opportunity to extend the existing n = 63 log /mpaq 
database to include more log 7mpaq values greater than 0.0. The log 7mhaq values for n = 18 phenols 
and their corresponding physicochemical properties that are necessary to determine their fit to RS were 



Pharmaceutics 2013, 5 



436 



published by Roberts, et al. in 1977 [10]. The fit of the n = 18 phenols to RS in the fit of the n = 62 
edited Flynn database to RS was published in 2007 [7]. In the edited Flynn database only n = 16 of the 
n = 62 molecules exhibited log 7mhaq values greater than 0.0 and of those n = 16, n = 11 were from 
among the n = 18 phenol subset [10]. Thus, phenols as a subset represent molecules that exhibit 
physicochemical properties (input) that give higher flux (output) than other types of molecules give. 

At present, only n = 6 of the n = 18 simple phenol subset from Roberts, et al. [10] have been 
included in the n = 63 log 7mpaq database and only n = 2 exhibit log 7mpaq greater than 0.0. In order to 
improve the correlation of the n = 63 log 7mpaq database with a matched n = 48 log 7mhaq database, 
the number of log 7mpaq and log 7mhaq greater than 0.0 in each database should be increased. Hence, 
n = 1 additional phenols have been selected from the n = 18 subset which exhibit physicochemical 
properties (input) for which RS predicts high log 7mpaq values (output). In addition, the n = 7 phenols 
exhibit an average log /mhaq value significantly greater than that of the n = 48 log /mhaq database: 
means + 95% confidence intervals of 0.04 + 0.42 log units and -1.06 + 0.31 log units, respectively. 
Given the increased range and total number of entries resulting from the addition of these n = 1 
compounds to the n = 48 log 7mhaq database and the n = 63 log /mpaq database, the fit of these 
databases to the RS should improve, and correlation of the log 7mpaq with log /mhaq values matched in 
these databases should also improve. 

Further, since the addition of the n = 1 new entries to the new n = 63 log 7mpaq database, each 
potentially exhibiting higher log 7mpaq values than the average of the initial n = 63 log 7mpaq values, 
will change the relative distribution of flux values in the database, it is imperative to determine if other 
models would then fit the database better than they did before the addition of the n = 1 new entries. 
Thus, we will also determine the fit of the new databases to the Kasting-Smith-Cooper (KSC) model [11] 
and to the Magnusson-Anissimov-Cross-Roberts (MACR) model [12] and compare these fits to the 
fit of RS to the new databases. 

2. Materials and Methods 

The phenolic compounds used are listed in Table 1. These compounds were obtained from Aldrich 
and their solubility values were acquired or approximated from literature sources. The phenols were all 
solids except for 3-methylphenol. 

The measurement of maximum flux through silicone was performed according to a literature 
procedure [13] at 32 °C, except that the silicone membrane was in contact with the receptor for only 
24 h to condition them. The receptor was a 7.1 pH phosphate buffer. 

The donor suspensions were prepared by stirring approximately 0.5 g (1 g in the case of 
3-methylphenol) of the compounds in 10 mL of water for 24 h. For all compounds, this surpassed the 
aqueous solubility by a factor of at least 20, which ensured saturation and excess solid/oil present in 
the donor phase. After the membranes were conditioned, the receptor phases were changed and the 
donor suspensions (first application, 1 mL) were applied; n = 3. The donor cells were sealed by 
Parafilm. Samples were taken from the receptor every 2-3 h after application. Following sample 
collection, the receptor phases were changed to ensure sink conditions, and the donor suspensions were 
either changed or had more solid/oil added to the existing suspension, depending upon the visible 
extent of depletion. After 4-5 sampling intervals, the donor suspensions were removed with methanol 
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and the receptors were changed. The membranes were leached with methanol in the donor phase for 
48-72 h with samples taken and receptor phases changed every 12-24 h to remove any residual phenol 
in the membrane. 

To ensure that flux data was not altered by possible membrane damage, a standard solute/solvent 
was applied and its flux determined. A donor suspension was prepared from 400 mg of theophylline 
suspended with stirring in 6 mL of propylene glycol (PG) for 24 h. This suspension (second 
application, 0.50 mL) was applied to all the silicone membranes after they were leached with methanol 
(see above). Samples were taken from the receptor every 24 h after application for at least 72 h so that 
at least 3 samples were obtained. Following sample collection, the receptor phases were changed and 
the donor suspensions were changed every other sampling interval. After 3-5 sampling intervals, 
the diffusion cells were disassembled and the membranes were placed in a methanol bath for 
maintenance leaching. 

The flux values of the first and second application were determined by UV absorption. The 
wavelengths (X E ) and molar absorptivities (s) used for the phenolic compounds are listed in Table 1. 
The log flux of theophylline through silicone from PG, log /mppg, for each membrane was found to be 
within the standard deviation of the literature value of -2.68 + 0.12 log units [13]. 

Nonlinear regression was performed by SPSS 20.0 (Rel. 20.0.0). The compounds were fitted to the 
Roberts-Sloan equation for maximum flux, log 7maq: 

log/MAQ = x + y log Sqct + (1 - y) log S AQ -zMW (2) 

to the KSC Equation: 

log/MAQ = x + y log Sqct _ z M w (3) 

and to the MACR Equation: 

log/MAQ = x — z MW (4) 

3. Results and Discussion 

The results are displayed in Table 1. All but 4-chloro-3,5-dimethylphenol exhibited a log 7mpaq 
greater than 0.0, and even it was very close. As a subset the n = 7 simple phenols had an average 
log /mpaq significantly greater than the average log 7mpaq of the n = 63 log /mpaq database: means + 
95% confidence intervals 1.03 + 0.45 log units and -0.42 + 0.29 log units, respectively. The average 
log 7mpaq in the n = 70 log /mpaq database has not significantly increased, but is no longer 
significantly less than 0.0: mean + 95% confidence interval -0.27 + 0.29. Unfortunately, the addition 
of the n = 7 phenols did not significantly increase the average log 7mhaq of the n = 55 log /mhaq 
database relative to the n = 48 log 7mhaq database: means + 95% confidence intervals, -0.918 ± 0.29 
log units and -1.058 + 0.31 log units, respectively. 
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Table 1. The relevant measured or literature physicochemical properties for the n = 7 
phenolic compounds used in this study. 



Cmpd. a 


MW 


Log S A q m 


Log ^OCT:AQ b 


Log Soct b ' d 






Log /mpaq c ' f 


Log /mhaq b ' f 


1 


143 


1.55 


3.10 


4.65 


283 


1241 


1.01 


0.29 


2 


157 


0.28 


3.39 


3.67 


285 


1041 


-0.027 


-0.95 


3 


122 


1.61 


2.35 


3.96 


277 


1668 


1.37 


0.17 


4 


108 


2.29 


1.95 


4.24 


276 


1614 


1.62 


0.53 


5 


163 


1.49 


3.08 


4.57 


285 


1791 


1.16 


0.27 


6 


197 


0.66 


3.69 


4.35 


312 


4518 


0.49 


-0.57 


7 


108 


2.36 


1.96 


4.32 


271 


1468 


1.61 


0.54 



a Substituted phenols. 1, 4-chloro-3-methyl; 2, 4-chloro-3,5-dimethyl; 3, 3,4-dimethyl; 4, 4-methyl; 5, 2,4-dichloro; 
6, 2,4,6-trichloro; 7, 3-methyl; b From Roberts et al. 1977 [10] and Majumdar et al. 2007 [7]; c Measured directly. 
d Solubility in water (S AQ ) or octanol (S 0 ct) m nmole cm -3 ; e Molar absorptivity coefficient in L mole -1 cm -1 ; f Maximum 
flux through silicone (7mpaq) or human stratum corneum (Jmhaq) from water in umole cm -2 If 1 . 



The addition of these n = 7 phenols to the n = 63 log /mpaq database and the n = 48 log 7mhaq 
database improved the fit of these databases to the RS as expected. The fit of the new n = 70 log 7mpaq 
database gave an r of 0.907, an average absolute residual log 7mpaq (Alog 7mpaq) of 0.300 log units 
and the coefficients x = -1.606, y = 0.695 and z = 0.00490 were all significant (p < 0.05): 

log/MPAQ = -1.606 + 0.695 logger + °- 305 lo g 5 AQ - 0.00490 MW (5) 

The fit of the n = 70 log /mpaq database is an improvement over the n = 63 log 7mpaq database, 
which had r = 0.896 and Alog 7mpaq = 0.310 log units, but had similar coefficient values: x = -1.607, 
y = 0.701, z = 0.00492. The fit of the new n = 55 log 7mhaq database gave an r 2 of 0.883, an average 
absolute residual log /mhaq (Alog 7mhaq) of 0.282 log units and the coefficients x = -3.005 and 
y = 0.654 were significant (p < 0.05), but the coefficient z = 0.001 12 was not significant (p = 0.25): 

log/MHAQ = -3.005 + 0.654 logger + 0346 lo g 5 AQ - 0.00112 MW (6) 

The lack of statistical significance for the z coefficient indicates a need to further extend the n = 55 
log 7mhaq database, since the significance of MW to maximum flux is well-established [12]. The fit of 
the n = 55 log 7mhaq database to the RS is an improvement over the n = 48 log /mhaq database, which 
had r = 0.867 and Alog 7mhaq = 0.331 log units and the coefficients x = -2.763, y = 0.635 and 
z = 0.00207. The x and y coefficients for the new log 7mhaq database are substantially closer to those 
coefficients determined for the n = 62 edited Flynn log 7mhaq database: x = -3.008, y = 0.732, 
z = 0.0048. Relevant results are displayed in Table 2 and Figures 1 and 2. Figure 1 shows a plot of 
experimental (Exp.) log 7mpaq versus log 7mpaq calculated (Calc.) from the coefficients for the fit of 
the n = 70 database to RS, and Figure 2 shows a plot of Exp. log /mhaq versus log 7mhaq Calc. from 
the coefficients for the fit of the n = 55 database to RS. 
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Table 2. The calculated (Calc.), predicted (Pred.), and experimental (Exp.) maximum flux 
values through silicone from water (log /mpaq) and through human stratum corneum from 
water (log /mhaq) for the n = 1 phenolic compounds. 



Cmpd. a 


Exp. Pred. n = 63 
log /mpaq b log /mpaq b ' c 


Calc. n = 70 

, j b,d 

'Og JMPAQ 


Exp. 
log /mhaq b 


Pred. n = 48 

log /mhaq b ' e 


Calc. n = 55 

log /mhaq b ' f 


1 


1.01 1.41 


1.40 


0.29 


0.46 


0.41 


2 


-0.027 0.28 


0.26 


-0.95 


-0.66 


-0.68 


3 


1.37 1.05 


1.04 


0.17 


0.087 


0.0053 


4 


1.62 1.52 


1.51 


0.53 


0.54 


0.44 


5 


1.16 1.24 


1.23 


0.27 


0.35 


0.32 


6 


0.49 0.67 


0.65 


-0.57 


-0.17 


-0.15 


7 


1.61 1.60 


1.59 


0.54 


0.62 


0.52 


Alog J MM} s 


0.200 


0.195 




0.159 


0.162 



a Substituted phenols. 1, 4-chloro-3-methyl; 2, 4-chloro-3,5-dimethyl; 3, 3,4-dimethyl; 4, 4-methyl; 5, 2,4-dichloro; 
6, 2,4,6-trichloro; 7, 3-methyl; b Given in units umole cm 2 h '; ' Using RS coefficients x = -1.607, y = 0.701, 
z = 0.00492 [6]; d Using RS coefficients x = -1.606, y = 0.695, z = 0.00490; e Using RS coefficients x = -2.763, y = 0.635 
and z = 0.00207 [6]; f Using RS coefficients x = -3.005, y = 0.654, z = 0.001 12; g Average absolute residual log 7 M paq or 
log 7 M haq f or the w = 7 phenols. 



Figure 1. The correlation of the calculated (Calc.) log 7mpaq from the fit of n = 70 to RS 
with the experimental (Exp.) log 7mpaq- The dashed lines represent the boundaries for residual 
log 7mpaq greater than 1.0, the solid line indicates points where the Calc. log 7mpaq is 
equivalent to the Exp. log 7mpaq- The filled circles indicate the n = 7 phenols. The Calc. 
log 7mpaq values were determined with Equation 5: log 7mpaq = -1.606 + 0.695 log 5oct + 
0.305 log Saq - 0.00490MW. r 2 = 0.907, average absolute residual log 7 MPA q = 0.300. 
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Figure 2. The correlation of the calculated (Calc.) log 7mhaq from the fit of n = 55 to RS 
with the experimental (Exp.) log /mhaq- The dashed lines represent the boundaries for 
residual Exp. log /mhaq greater than 1.0, and the solid line indicates points where the 
Calc. log /mhaq is equivalent to the Exp. log /mhaq- The filled circles indicate the n = 1 
phenols. The Calc. log 7mhaq values were determined with Equation 6: log 7mhaq = -3.005 
+ 0.654 log Soct + 0.346 log S A q ~ 0.00112 MW, r 2 = 0.883, average absolute residual 
log Jmhaq = 0.282. 
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Plots of the individual independent variables, log 5oct, log 5aq and MW, against log 7mpaq in = 70) 

2 2 

and against log 7mhaq in = 55) gave the following r values: (a) log 7mpaq versus log Soct, r = 0.677; 

2 2 2 

versus log 5aq, r = 0.554; versus MW, r = 0.541; (b) log 7mhaq versus log Soct, r = 0.603; versus 

2 2 

log 5aq, r = 0.526; versus MW, r = 0.520. All regression equations had statistically significant 
ip < 0.05) slope and intercept estimates. In each case, the best regression of the individual independent 
variables against flux values was by log 5oct- It is worth noting that the n = 1 phenol subset gives the 
following r values and significance profiles when the individual independent variables log 5oct and 
log 5aq are plotted against log /mpaq and log /mhaq: (a) log /mpaq versus log Soct, r = 0.180, without 
a statistically significant slope ip = 0.34) or intercept ip = 0.51); versus log 5 A q, r = 0.949, with a 
statistically significant slope ip < 0.05), but without a statistically significant intercept ip = 0.50); 
(b) log /mhaq versus log Soct, r 2 = 0.335, without a statistically significant slope ip = 0.17) or intercept 
(p = 0.18); versus log Saq, r 2 = 0.933, with statistically significant ip < 0.05) slope and intercept. The 
improved dependence of maximum flux from water on the aqueous solubility of highly water-soluble 
compounds is demonstrated here as a property of both silicone and human stratum corneum, and will 
be a topic of future investigations. When two individual independent variables, log Soct and MW, 
from the n = 10 log 7mpaq database were fitted to the KSC equation (Equation 3) the following x, y and 
z coefficients to the parameters were obtained along with r 2 and Alog /mpaq values: x = -0.923, 
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y = 0.794, z = 0.0089, r = 0.797 and Alog 7 M paq = 0.431. All but the y coefficient (p = 0.069) were 
statistically significant (p < 0.05): 



log/MPAQ = -0.923 + 0.794 logS OCT - 0.0089 MW 



(7) 



Figure 3 shows a plot of Exp. log /mpaq versus log /mpaq values calculated from the coefficients for 
the fit of the n = 70 database to Equation 7. Although the r was substantially improved by including 

2 2 

MW with log Soct as independent variables in the regression against flux, the r was poorer than the r 
for the fit of all three independent variables to the Roberts-Sloan Equation (Equation 5). Similarly, 
when two individual independent variables, log Soct and MW, from the n = 55 log 7mhaq database 
were fit to the KSC equation (Equation 8) the following x, y and z coefficients to the parameters were 

2 2 

obtained along with r and the Alog 7mhaq values: x = -1.252, y = 0.602, z = 0.0080, r = 0.723 and 
Alog 7mhaq = 0.441. The estimates for the coefficients were all statistically significant (p < 0.05): 

log/ MHAQ = -1.252 + 0.602 logS OCT - 0.0080 MW (8) 

Figure 3. The correlation of the calculated (Calc.) log 7mpaq from the fit of n = 70 to KSC 
with the experimental (Exp.) log 7mpaq- The dashed lines represent the boundaries for 
residual Exp. log 7mpaq greater than 1.0, the solid line indicates points where the Calc. 
log 7mpaq is equivalent to the Exp. log 7mpaq- The filled circles indicate the n = 1 phenols. 
The Calc. log 7mpaq values were determined with Equation 7: log/ MPA Q = —0.923 + 

0.794 logSocT ~ °- 0089 MW > r 2 = 0-797, average absolute residual log 7 M paq = 0.431. 
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Figure 4 shows a plot of Exp. log 7mhaq versus log 7mhaq values calculated from the coefficients 
for the fit of the n = 55 database to Equation 8. Again, although a substantial improvement in r 1 was 
obtained by including MW with log 5oct as independent variables in the regression against flux, the r 
was poorer than the r for the fit of all three independent variables to the RS equation (Equation 6). 
The fit of both databases to the MACR equation (Equation 4), which is the remaining model used to 
predict maximum flux, is simply the regression of MW against log 7mpaq or log /mhaq shown above to 
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give a somewhat poorer fit than regression of the two individual independent variables, log Soct or log Saq, 
against log 7mpaq or log 7mhaq- Finally, it should be noted the popular Potts-Guy Equation [14] was 
not included as a model because its output is permeability coefficient which is not clinically relevant. 

Figure 4. The correlation of the calculated (Calc.) log 7mhaq from the fit of n = 55 to KSC 
with the experimental (Exp.) log 7mhaq- The dashed lines represent the boundaries for 
residual Exp. log /mhaq greater than 1.0, and the solid line indicates points where the Calc. 
log 7mhaq is equivalent to the Exp. log 7mhaq- The filled circles indicate the n = 7 phenols. 
The Calc. log 7mhaq values were determined with Equation 8: log/ MH AQ = —1-252 + 
0.602 logSocT ~ °- 0080 MW > r 2 = 0-723, average absolute residual log 7 M haq = 0.441. 
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The new n = 52 database of compounds contributing to both the n = 70 log 7mpaq database and the 
n = 55 log 7mhaq database also gives a higher correlation between Exp. log /mpaq and Exp. log 7mhaq 
than the previous n = 45 database. A linear regression yielded the expression Exp. log 7mhaq = 0.859 
Exp. log 7mpaq ~ 0.837, r 2 = 0.856, which is an improvement over r 2 = 0.838 for the n = 45 database. 
Figure 5 shows the plot of Exp. log 7mpaq versus Exp. log 7mhaq for n = 52. 
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Figure 5. The correlation of the n = 52 log 7mhaq with log 7mpaq database. The dashed 
lines represent the boundaries for residual Exp. log /mhaq greater than 1.0, and the solid 
line is the regression equation. The filled circles indicate the n = 7 phenols. The regression 
information is in the figure. 
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4. Conclusions 

The addition of the n = 1 phenols improved all aspects of the log 7mpaq database with regards to 
comparison with the matched log 7mhaq database. Along with strengthening the validity of the 
assertion that silicone membranes are good surrogates for human stratum corneum, these improvements 
demonstrate that this surrogate nature holds for a wider range of log 7mpaq and log 7mhaq values than 
had been previously reported. 
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